function [varargout]=Equ_NK(PP,SS,order)

%'Eta_Z','Eta_Zp','Eta_Z_N','Eta_Z_Np','Eta_Z_H','Eta_Z_Hp','Eta_M_dom','Eta_M_domp',...
 %'Eta_M_ext','Eta_M_extp','Eta_Y_H','Eta_Y_Hp','N_N','N_H','w','mc_N',...
 %'mc_H','Y_N','Y_H','Profit_N','Profit_H','PAC_N','PAC_H','Y_Np',...
 %'Y_Hp','L','C_N','C_T','C_H','C_F','p_T','p_N',...
 %'p_H','p_F','Pir','Pirp','Pir_N','Pir_H','Pir_F','Lag_p_T',...
 %'Lag_p_N','Lag_p_H','Lag_p_F','Lag_p_Tp','Lag_p_Np','Lag_p_Hp','Lag_p_Fp','p_Np',...
 %'p_Hp','Pir_Np','Pir_Hp','ir_dom','ir_ext','rr_dom','rr_ext','ir',...
 %'TAU','B','T','ir_s','Pir_H_s','Pir_F_s','dE','dEp',...
 %'Y','C','Profit','PAC','GB','Lag_GB','Lag_GBp','Lag_B',...
 %'Lag_Bp','Lag_ir','Lag_irp',
 
 PP_RHO_Z=PP.RHO_Z;
PP_RHO_Z_N=PP.RHO_Z_N;
PP_RHO_Z_H=PP.RHO_Z_H;
PP_RHO_M_dom=PP.RHO_M_dom;
PP_RHO_M_ext=PP.RHO_M_ext;
PP_RHO_Y_H=PP.RHO_Y_H;
PP_ALPHA_N=PP.ALPHA_N;
PP_ALPHA_H=PP.ALPHA_H;
PP_EPSILON_N=PP.EPSILON_N;
PP_EPSILON_H=PP.EPSILON_H;
PP_THETA_N=PP.THETA_N;
PP_THETA_H=PP.THETA_H;
PP_Cons_FracT=PP.Cons_FracT;
PP_Cons_FracH=PP.Cons_FracH;
PP_Cons_ElasTN=PP.Cons_ElasTN;
PP_Cons_ElasHF=PP.Cons_ElasHF;
PP_Taylor_Pi=PP.Taylor_Pi;
PP_Taylor_Pi_H=PP.Taylor_Pi_H;
PP_Taylor_Pi_N=PP.Taylor_Pi_N;
PP_Taylor_dE=PP.Taylor_dE;
PP_Taylor_ir=PP.Taylor_ir;
PP_Flag_dE=PP.Flag_dE;
PP_Flag_Price=PP.Flag_Price;

 SS_Z_N=SS.Z_N;
SS_Z_H=SS.Z_H;
SS_N_N=SS.N_N;
SS_N_H=SS.N_H;
SS_w=SS.w;
SS_mc_N=SS.mc_N;
SS_mc_H=SS.mc_H;
SS_Y_N=SS.Y_N;
SS_Y_H=SS.Y_H;
SS_Profit_N=SS.Profit_N;
SS_Profit_H=SS.Profit_H;
SS_PAC_N=SS.PAC_N;
SS_PAC_H=SS.PAC_H;
SS_C_T=SS.C_T;
SS_C_N=SS.C_N;
SS_C_H=SS.C_H;
SS_C_F=SS.C_F;
SS_L=SS.L;
SS_p_T=SS.p_T;
SS_p_N=SS.p_N;
SS_p_H=SS.p_H;
SS_p_F=SS.p_F;
SS_Pir=SS.Pir;
SS_Pir_N=SS.Pir_N;
SS_Pir_H=SS.Pir_H;
SS_Pir_F=SS.Pir_F;
SS_ir_dom=SS.ir_dom;
SS_ir_ext=SS.ir_ext;
SS_ir=SS.ir;
SS_TAU=SS.TAU;
SS_B=SS.B;
SS_T=SS.T;
SS_Y=SS.Y;
SS_C=SS.C;
SS_Profit=SS.Profit;
SS_PAC=SS.PAC;
SS_G=SS.G;
SS_GB=SS.GB;
SS_C_H_s=SS.C_H_s;

 Eta_Z=0;
Eta_Zp=0;
Eta_Z_N=0;
Eta_Z_Np=0;
Eta_Z_H=0;
Eta_Z_Hp=0;
Eta_M_dom=0;
Eta_M_domp=0;
Eta_M_ext=0;
Eta_M_extp=0;
Eta_Y_H=0;
Eta_Y_Hp=0;
N_N=0;
N_H=0;
w=0;
mc_N=0;
mc_H=0;
Y_N=0;
Y_H=0;
Profit_N=0;
Profit_H=0;
PAC_N=0;
PAC_H=0;
Y_Np=0;
Y_Hp=0;
L=0;
C_N=0;
C_T=0;
C_H=0;
C_F=0;
p_T=0;
p_N=0;
p_H=0;
p_F=0;
Pir=0;
Pirp=0;
Pir_N=0;
Pir_H=0;
Pir_F=0;
Lag_p_T=0;
Lag_p_N=0;
Lag_p_H=0;
Lag_p_F=0;
Lag_p_Tp=0;
Lag_p_Np=0;
Lag_p_Hp=0;
Lag_p_Fp=0;
p_Np=0;
p_Hp=0;
Pir_Np=0;
Pir_Hp=0;
ir_dom=0;
ir_ext=0;
rr_dom=0;
rr_ext=0;
ir=0;
TAU=0;
B=0;
T=0;
ir_s=0;
Pir_H_s=0;
Pir_F_s=0;
dE=0;
dEp=0;
Y=0;
C=0;
Profit=0;
PAC=0;
GB=0;
Lag_GB=0;
Lag_GBp=0;
Lag_B=0;
Lag_Bp=0;
Lag_ir=0;
Lag_irp=0;
if order==0
res=ones(51,1);
res(1)=Eta_Z*PP_RHO_Z - Eta_Zp;
res(2)=Eta_Z_N*PP_RHO_Z_N - Eta_Z_Np;
res(3)=Eta_Z_H*PP_RHO_Z_H - Eta_Z_Hp;
res(4)=Eta_M_dom*PP_RHO_M_dom - Eta_M_domp;
res(5)=Eta_M_ext*PP_RHO_M_ext - Eta_M_extp;
res(6)=Eta_Y_H*PP_RHO_Y_H - Eta_Y_Hp;
res(7)=SS_Z_N*exp(Eta_Z + Eta_Z_N)*(SS_N_N*exp(N_N))^(1 - PP_ALPHA_N) - SS_Y_N*exp(Y_N);
res(8)=SS_Z_H*exp(Eta_Z + Eta_Z_H)*(SS_N_H*exp(N_H))^(1 - PP_ALPHA_H) - SS_Y_H*exp(Y_H);
res(9)=- SS_mc_N*exp(mc_N) - (SS_w*SS_N_N*exp(-Y_N)*exp(N_N + w))/(SS_Y_N*(PP_ALPHA_N - 1));
res(10)=- SS_mc_H*exp(mc_H) - (SS_w*SS_N_H*exp(-Y_H)*exp(N_H + w))/(SS_Y_H*(PP_ALPHA_H - 1));
res(11)=(Pir_Np*(SS_Pir_N + 1))/(SS_ir + 1) - Pir_N + (PP_EPSILON_N*(mc_N - p_N))/(PP_THETA_N*(SS_Pir_N + 1));
res(12)=((Pir_Hp - PP_Flag_Price*dEp*(SS_Pir_H + 1))*(SS_Pir_H + 1))/(SS_ir + 1) - Pir_H + PP_Flag_Price*dE*(SS_Pir_H + 1) + (PP_EPSILON_H*(mc_H - p_H))/(PP_THETA_H*(SS_Pir_H + 1));
res(13)=-PAC_N;
res(14)=-PAC_H;
res(15)=SS_Y_N*SS_p_N*exp(Y_N + p_N) - Profit_N - SS_PAC_N - SS_Profit_N - SS_w*SS_N_N*exp(N_N + w) - PAC_N;
res(16)=SS_Y_H*SS_p_H*exp(Y_H + p_H) - Profit_H - SS_PAC_H - SS_Profit_H - SS_w*SS_N_H*exp(N_H + w) - PAC_H;
res(17)=(PP_Cons_FracT*SS_C*exp(C))/(SS_p_T*exp(p_T))^PP_Cons_ElasTN - SS_C_T*exp(C_T);
res(18)=- SS_C_N*exp(C_N) - (SS_C*exp(C)*(PP_Cons_FracT - 1))/(SS_p_N*exp(p_N))^PP_Cons_ElasTN;
res(19)=(PP_Cons_FracH*SS_C_T*exp(C_T))/((SS_p_H*exp(p_H - p_T))/SS_p_T)^PP_Cons_ElasHF - SS_C_H*exp(C_H);
res(20)=- SS_C_F*exp(C_F) - (SS_C_T*exp(C_T)*(PP_Cons_FracH - 1))/((SS_p_F*exp(p_F - p_T))/SS_p_T)^PP_Cons_ElasHF;
res(21)=PP_Cons_FracT*p_T - p_N*(PP_Cons_FracT - 1);
res(22)=PP_Cons_FracH*p_H - p_T - p_F*(PP_Cons_FracH - 1);
res(23)=exp(p_N - Lag_p_N)*(Pir + SS_Pir + 1) - SS_Pir_N - Pir_N - 1;
res(24)=exp(p_H - Lag_p_H)*(Pir + SS_Pir + 1) - SS_Pir_H - Pir_H - 1;
res(25)=exp(p_F - Lag_p_F)*(Pir + SS_Pir + 1) - SS_Pir_F - Pir_F - 1;
res(26)=Eta_M_dom - ir + PP_Taylor_Pi*Pir + PP_Taylor_Pi_H*Pir_H + PP_Taylor_Pi_N*Pir_N + PP_Taylor_dE*dE + Lag_ir*PP_Taylor_ir;
res(27)=(GB + SS_GB)/(SS_ir + ir + 1) - SS_T - T - (Lag_GB + SS_GB)/(Pir + SS_Pir + 1) - SS_G + SS_L*SS_w*exp(L + w)*(SS_TAU + TAU);
res(28)=-GB;
res(29)=-T;
res(30)=ir - ir_dom;
res(31)=ir_s - ir_ext;
res(32)=ir_dom - Pirp - rr_dom;
res(33)=dEp - Pirp + ir_ext - rr_ext;
res(34)=Eta_M_ext - ir_s;
res(35)=SS_C_H*exp(C_H) - SS_Y_H*exp(Y_H) + exp(PP_Cons_ElasHF*(p_F - p_H))*(SS_C_H_s + Eta_Y_H*SS_Y_H);
res(36)=Pir_F_s - Pir_F + dE;
res(37)=-Pir_H_s;
res(38)=-Pir_F_s;
res(39)=dEp - ir_dom + ir_ext;
res(40)=SS_C_N*exp(C_N) - SS_Y_N*exp(Y_N);
res(41)=SS_N_H*exp(N_H) - SS_L*exp(L) + SS_N_N*exp(N_N);
res(42)=SS_Y_H*SS_p_H*exp(Y_H + p_H) - SS_Y*exp(Y) + SS_Y_N*SS_p_N*exp(Y_N + p_N);
res(43)=PAC_H - PAC + PAC_N - SS_PAC + SS_PAC_H + SS_PAC_N;
res(44)=Profit_H - Profit + Profit_N - SS_Profit + SS_Profit_H + SS_Profit_N;
res(45)=p_N - Lag_p_Np;
res(46)=p_T - Lag_p_Tp;
res(47)=p_H - Lag_p_Hp;
res(48)=p_F - Lag_p_Fp;
res(49)=GB - Lag_GBp;
res(50)=B - Lag_Bp;
res(51)=ir - Lag_irp;
varargout{1}=res;
elseif order==1
fir_ord.Eta_Z=[PP_RHO_Z;0;0;0;0;0;SS_Z_N*exp(Eta_Z + Eta_Z_N)*(SS_N_N*exp(N_N))^(1 - PP_ALPHA_N);SS_Z_H*exp(Eta_Z + Eta_Z_H)*(SS_N_H*exp(N_H))^(1 - PP_ALPHA_H);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Zp=[-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Z_N=[0;PP_RHO_Z_N;0;0;0;0;SS_Z_N*exp(Eta_Z + Eta_Z_N)*(SS_N_N*exp(N_N))^(1 - PP_ALPHA_N);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Z_Np=[0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Z_H=[0;0;PP_RHO_Z_H;0;0;0;0;SS_Z_H*exp(Eta_Z + Eta_Z_H)*(SS_N_H*exp(N_H))^(1 - PP_ALPHA_H);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Z_Hp=[0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_M_dom=[0;0;0;PP_RHO_M_dom;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_M_domp=[0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_M_ext=[0;0;0;0;PP_RHO_M_ext;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_M_extp=[0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Y_H=[0;0;0;0;0;PP_RHO_Y_H;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_Y_H*exp(PP_Cons_ElasHF*(p_F - p_H));0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Eta_Y_Hp=[0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.N_N=[0;0;0;0;0;0;-(SS_N_N*SS_Z_N*exp(Eta_Z + Eta_Z_N)*exp(N_N)*(PP_ALPHA_N - 1))/(SS_N_N*exp(N_N))^PP_ALPHA_N;0;-(SS_w*SS_N_N*exp(-Y_N)*exp(N_N + w))/(SS_Y_N*(PP_ALPHA_N - 1));0;0;0;0;0;-SS_w*SS_N_N*exp(N_N + w);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_N_N*exp(N_N);0;0;0;0;0;0;0;0;0;0];
fir_ord.N_H=[0;0;0;0;0;0;0;-(SS_N_H*SS_Z_H*exp(Eta_Z + Eta_Z_H)*exp(N_H)*(PP_ALPHA_H - 1))/(SS_N_H*exp(N_H))^PP_ALPHA_H;0;-(SS_w*SS_N_H*exp(-Y_H)*exp(N_H + w))/(SS_Y_H*(PP_ALPHA_H - 1));0;0;0;0;0;-SS_w*SS_N_H*exp(N_H + w);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_N_H*exp(N_H);0;0;0;0;0;0;0;0;0;0];
fir_ord.w=[0;0;0;0;0;0;0;0;-(SS_w*SS_N_N*exp(-Y_N)*exp(N_N + w))/(SS_Y_N*(PP_ALPHA_N - 1));-(SS_w*SS_N_H*exp(-Y_H)*exp(N_H + w))/(SS_Y_H*(PP_ALPHA_H - 1));0;0;0;0;-SS_w*SS_N_N*exp(N_N + w);-SS_w*SS_N_H*exp(N_H + w);0;0;0;0;0;0;0;0;0;0;SS_L*SS_w*exp(L + w)*(SS_TAU + TAU);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.mc_N=[0;0;0;0;0;0;0;0;-SS_mc_N*exp(mc_N);0;PP_EPSILON_N/(PP_THETA_N*(SS_Pir_N + 1));0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.mc_H=[0;0;0;0;0;0;0;0;0;-SS_mc_H*exp(mc_H);0;PP_EPSILON_H/(PP_THETA_H*(SS_Pir_H + 1));0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Y_N=[0;0;0;0;0;0;-SS_Y_N*exp(Y_N);0;(SS_w*SS_N_N*exp(-Y_N)*exp(N_N + w))/(SS_Y_N*(PP_ALPHA_N - 1));0;0;0;0;0;SS_Y_N*SS_p_N*exp(Y_N + p_N);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_Y_N*exp(Y_N);0;SS_Y_N*SS_p_N*exp(Y_N + p_N);0;0;0;0;0;0;0;0;0];
fir_ord.Y_H=[0;0;0;0;0;0;0;-SS_Y_H*exp(Y_H);0;(SS_w*SS_N_H*exp(-Y_H)*exp(N_H + w))/(SS_Y_H*(PP_ALPHA_H - 1));0;0;0;0;0;SS_Y_H*SS_p_H*exp(Y_H + p_H);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_Y_H*exp(Y_H);0;0;0;0;0;0;SS_Y_H*SS_p_H*exp(Y_H + p_H);0;0;0;0;0;0;0;0;0];
fir_ord.Profit_N=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0];
fir_ord.Profit_H=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0];
fir_ord.PAC_N=[0;0;0;0;0;0;0;0;0;0;0;0;-1;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0];
fir_ord.PAC_H=[0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0];
fir_ord.Y_Np=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Y_Hp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.L=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_L*SS_w*exp(L + w)*(SS_TAU + TAU);0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_L*exp(L);0;0;0;0;0;0;0;0;0;0];
fir_ord.C_N=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_C_N*exp(C_N);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_C_N*exp(C_N);0;0;0;0;0;0;0;0;0;0;0];
fir_ord.C_T=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_C_T*exp(C_T);0;(PP_Cons_FracH*SS_C_T*exp(C_T))/((SS_p_H*exp(p_H - p_T))/SS_p_T)^PP_Cons_ElasHF;-(SS_C_T*exp(C_T)*(PP_Cons_FracH - 1))/((SS_p_F*exp(p_F - p_T))/SS_p_T)^PP_Cons_ElasHF;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.C_H=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_C_H*exp(C_H);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_C_H*exp(C_H);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.C_F=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_C_F*exp(C_F);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.p_T=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-(PP_Cons_FracT*PP_Cons_ElasTN*SS_C*SS_p_T*exp(C)*exp(p_T))/(SS_p_T*exp(p_T))^(PP_Cons_ElasTN + 1);0;(PP_Cons_FracH*PP_Cons_ElasHF*SS_C_T*SS_p_H*exp(p_H - p_T)*exp(C_T))/(SS_p_T*((SS_p_H*exp(p_H - p_T))/SS_p_T)^(PP_Cons_ElasHF + 1));-(PP_Cons_ElasHF*SS_C_T*SS_p_F*exp(p_F - p_T)*exp(C_T)*(PP_Cons_FracH - 1))/(SS_p_T*((SS_p_F*exp(p_F - p_T))/SS_p_T)^(PP_Cons_ElasHF + 1));PP_Cons_FracT;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0];
fir_ord.p_N=[0;0;0;0;0;0;0;0;0;0;-PP_EPSILON_N/(PP_THETA_N*(SS_Pir_N + 1));0;0;0;SS_Y_N*SS_p_N*exp(Y_N + p_N);0;0;(PP_Cons_ElasTN*SS_C*SS_p_N*exp(C)*exp(p_N)*(PP_Cons_FracT - 1))/(SS_p_N*exp(p_N))^(PP_Cons_ElasTN + 1);0;0;1 - PP_Cons_FracT;0;exp(p_N - Lag_p_N)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_Y_N*SS_p_N*exp(Y_N + p_N);0;0;1;0;0;0;0;0;0];
fir_ord.p_H=[0;0;0;0;0;0;0;0;0;0;0;-PP_EPSILON_H/(PP_THETA_H*(SS_Pir_H + 1));0;0;0;SS_Y_H*SS_p_H*exp(Y_H + p_H);0;0;-(PP_Cons_FracH*PP_Cons_ElasHF*SS_C_T*SS_p_H*exp(p_H - p_T)*exp(C_T))/(SS_p_T*((SS_p_H*exp(p_H - p_T))/SS_p_T)^(PP_Cons_ElasHF + 1));0;0;PP_Cons_FracH;0;exp(p_H - Lag_p_H)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;-PP_Cons_ElasHF*exp(PP_Cons_ElasHF*(p_F - p_H))*(SS_C_H_s + Eta_Y_H*SS_Y_H);0;0;0;0;0;0;SS_Y_H*SS_p_H*exp(Y_H + p_H);0;0;0;0;1;0;0;0;0];
fir_ord.p_F=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;(PP_Cons_ElasHF*SS_C_T*SS_p_F*exp(p_F - p_T)*exp(C_T)*(PP_Cons_FracH - 1))/(SS_p_T*((SS_p_F*exp(p_F - p_T))/SS_p_T)^(PP_Cons_ElasHF + 1));0;1 - PP_Cons_FracH;0;0;exp(p_F - Lag_p_F)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;PP_Cons_ElasHF*exp(PP_Cons_ElasHF*(p_F - p_H))*(SS_C_H_s + Eta_Y_H*SS_Y_H);0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0];
fir_ord.Pir=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;exp(p_N - Lag_p_N);exp(p_H - Lag_p_H);exp(p_F - Lag_p_F);PP_Taylor_Pi;(Lag_GB + SS_GB)/(Pir + SS_Pir + 1)^2;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pirp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_N=[0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;PP_Taylor_Pi_N;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_H=[0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;-1;0;PP_Taylor_Pi_H;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_F=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_p_T=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_p_N=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-exp(p_N - Lag_p_N)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_p_H=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-exp(p_H - Lag_p_H)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_p_F=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-exp(p_F - Lag_p_F)*(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_p_Tp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0];
fir_ord.Lag_p_Np=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0];
fir_ord.Lag_p_Hp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0];
fir_ord.Lag_p_Fp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0];
fir_ord.p_Np=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.p_Hp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_Np=[0;0;0;0;0;0;0;0;0;0;(SS_Pir_N + 1)/(SS_ir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_Hp=[0;0;0;0;0;0;0;0;0;0;0;(SS_Pir_H + 1)/(SS_ir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.ir_dom=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;1;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.ir_ext=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;1;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.rr_dom=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.rr_ext=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.ir=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;-(GB + SS_GB)/(SS_ir + ir + 1)^2;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1];
fir_ord.TAU=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;SS_L*SS_w*exp(L + w);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.B=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0];
fir_ord.T=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.ir_s=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_H_s=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Pir_F_s=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;-1;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.dE=[0;0;0;0;0;0;0;0;0;0;0;PP_Flag_Price*(SS_Pir_H + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;PP_Taylor_dE;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.dEp=[0;0;0;0;0;0;0;0;0;0;0;-(PP_Flag_Price*(SS_Pir_H + 1)^2)/(SS_ir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0;0;0;0;1;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Y=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-SS_Y*exp(Y);0;0;0;0;0;0;0;0;0];
fir_ord.C=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;(PP_Cons_FracT*SS_C*exp(C))/(SS_p_T*exp(p_T))^PP_Cons_ElasTN;-(SS_C*exp(C)*(PP_Cons_FracT - 1))/(SS_p_N*exp(p_N))^PP_Cons_ElasTN;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Profit=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0];
fir_ord.PAC=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0;0;0;0;0;0;0];
fir_ord.GB=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1/(SS_ir + ir + 1);-1;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;1;0;0];
fir_ord.Lag_GB=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1/(Pir + SS_Pir + 1);0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_GBp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0;0];
fir_ord.Lag_B=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_Bp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1;0];
fir_ord.Lag_ir=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;PP_Taylor_ir;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0];
fir_ord.Lag_irp=[0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;-1];
varargout{1}=fir_ord;

 end
